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ABSTRACT 



We discuss the y-ray signal to be expected from dark matter (DM) annihilations at the Galactic Center. To describe the DM distribution 
in the Galactic halo we base on the Jeans equation for self-gravitating, anisotropic equilibria. In solving the Jeans equation, we adopt 
the specific correlation between the density p(r) and the velocity dispersion cr 2 (r) expressed by the powerlaw behavior of the DM 
'entropy' K = cr 2 /p 2/3 oc r" with a ss 1.25 - 1.3. Indicated (among others) by several recent A'-body simulations, this correlation 
is privileged by the form of the radial pressure term in the Jeans equation, and yields a main body profile consistent with the classic 
self-similar development of DM halos. In addition, we require the Jeans solutions to satisfy regular boundary conditions both at the 
center (finite pressure, round gravitational potential) and in the outskirts (finite overall mass). With these building blocks we derive 
physical solutions, dubbed '^-profiles'. We find the one with a = 1.25, suitable for the Galaxy halo, to be intrinsically flatter at 
the center relative to the empirical NFW formula, yet steeper than the empirical Einasto profile. So on scales of 10 _1 deg it yields 
annihilation fluxes lower by a factor 5 than the former yet higher by a factor 10 than the latter; such fluxes will eventually fall within 
the reach of the Fermi satellite. We show the effectiveness of the a-profile in relieving the astrophysical uncertainties related to the 
macroscopic DM distribution, and discuss its expected performance as a tool instrumental to interpret the upcoming y-ray data in 
terms of DM annihilation. 

Key words, dark matter - gamma rays: observations - galaxies: evolution - Galaxy: halo - methods: analytical 



1. Introduction 

Several astrophysic al and cosmological probes (for a review see 
IBertone et al.l20 05) have firmly established that baryons - which 



stars, planets, and (known) living creatures are made of - consti- 
tute only some 15% of the total matter content in the Universe 
adding to the dominant dark energy component. The rest is in the 
form of 'cold dark matter' (DM), i.e., massive particles that were 
non-relativistic at decoupling, do not emit/absorb radiations, and 
basically do not interact with themselves and with the baryons 
except via long-range gravitational forces. 

However, no 'direct' detection o f the DM has been made so 
far, other than [Bernabei et all (120081) . Thus the microscopic na- 
ture of the DM remains largely a mystery; several clues suggest 
as a promising candidate or component the lig htest supersym - 
metric particle, the 'neutralino' (for a review see Bertone 2009). 
Given that the latter' s mass, depending on the specific super- 
symmetric model, ranges from several GeVs to tens of TeVs, its 
laboratory production requires an accelerator at le ast as power- 
ful as the newly-born Large Hadron Collider (see Baer & Tata 
2009); the discovery of supersymmetry and specifically of the 
neutralino is one of the main aims for the current experiments in 
high-energy physics. 

Meanwhile, evidence for the DM can be looked for 'indi- 
rectly' in the sky. In fact, the basic aims of the recently launched 
Fermi satellite include the search for y-ray signals due to the 
annihilation of DM particles at the Galactic Center (GC) and in 
nearby galaxies (see discussion in § 4). The former provides a fa- 
vorable target being closest to us, with the DM density expected 
to increase in moving toward the inner regions of a galaxy. 



However, the GC is also a crowded region, and it remains a chal- 
lenging task to separate the DM signal from the contributions 
of other astrophysical sources and backgrounds whose energy 
spectrum and angular distribution are poorly known. 

In principle, if one can predict the strength and angular dis- 
tribution of the annihilation signal itself, then the y-ray observa- 
tions would elicit, or put 'indirect' constraints on the (combined) 
microscopic properties of the DM particles like mass, annihila- 
tion cross section and channels. This approach has been pursued 
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extensively (e.g., Bergstrom et al. 1998; Fornengo et al. 2004; 
IStrigaril l2007; Bert one et alj|2009tlsen>ico & HoopeJl2009l) bu t 
suffers yet of large uncertainties (see ICesarini et alj 12004). 
mainly related to the poor knowledge of the macroscopic DM 
distribution p(r) throughout the Galaxy. 

Since the annihilation rate scales like p 2 (r), such uncertain- 
ties are maximized near the center right where detection is fa- 
vored. Note that similar if milder uncertainties affect the source 
function of the electrons originated from DM annihilations by 
production or cascading; these diffuse outwards and interact 
with the Galactic magnetic field and with the interstellar light 
to produce synchrot ron emission observed in the radio band (see 
IBertone et all 2009). and inverse Compton radiation observable 



tsertone et al. ZUUV ), ana inverse Comptc 
in y rays (see Papucci & Strumia 2009). 



Traditionally, the density profile of an equilibrium DM struc- 
ture, or 'halo', is rendered in terms of different empirical formu- 
las that fit the results of TV-body simulations and to some ex- 
tent the stellar observations. Perhaps the most popular one is 
the N avarro, Frenk & White (hereafter NFW; see iNavarro et al.l 
1997) profile, that has an asymptotic inner slope p(r) oc r , 
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goes over to a powerlaw behavior p(r) oc r~ 2 in the halo's mid- 
dle, and declines as p(r) oc r~ 3 in the outer regions. Despite its 
widespread use in the literature, clearly this expression cannot 
account for the actual DM distribution in the inner regions of a 
galaxy halo where it would imply a centrally angled gravitational 
potential well and an infinite pressure, nor in the halo outskirts 
where it would yield a diverging overall mass. 

Other empirical density profiles have been proposed but suf- 
fer of similarly unph ysical features; e.g., the Moore profile (see 
iDiemand et al.ll2005l) goes like p(r) oc r -1 ' 2 and implies a grav- 
itational force diverging towards the center, while the Einasto 
profile (see Graham et al. 2006) behaves like p(r) oc e~ , so it 
yields a vanishing pressure there. We stress that the differences 
in the predicted annihilation signals under these DM distribu- 
tions turn out to be quite considerable; for example, the ratio of 
the NFW to the Einasto squared density averaged over 1 degree 
(about 150 pc) comes to a factor 10 when normalized at the Sun's 
location (see also discussion in § 4). 

Our stand here is that the macroscopic uncertainties yield- 
ing such differences can, and ought to be relieved. To this pur- 
pose, in § 2 we present the physical density distributions that 
we dub a-profiles; these are solutions of the Jeans equation that 
satisfy regular inner and outer boundary conditions. In § 3 we 
use the a-profile suitable for the Galaxy halo as the macroscopic 
benchmark to evaluate the DM annihilation signal expected from 
the GC. As for the microscopic sector, we base on a standard 
model for the mass, cross section and annihilation channel of 
the DM particles, the extension to more complex microphysics 
being straightforward. Finally, our findings are summarized and 
discussed in § 4. 

Throughout this wor k we adopt a standard, flat cosmology 
(see|Dunklev et al. 2009) with normalized matter density Q.m — 
0.27, and Hubble constant H = 72 km s -1 Mpc~'. 



2. Development and structure of DM halos 

Galaxies are widely held to form under the drive of the gravita- 
tional instability that acts on initial perturbations modulating the 
cosmic density of the dominant cold DM component. At first the 
instability is kept in check by the cosmic expansion, but when 
the local gravity prevails collapse sets in, and form a DM halo 
in equilibrium under self-gravity. The amplitude of more mas- 
sive perturbations is smaller, so the formation is progressive in 
time and hierarchic al in mass, with the largest structures forming 
typically later (see Peeb les et alJI 1983L for a review). 



2. 1 . Two-stage evolution 

Such a formation history has been resolved t o a consider- 
able detail by man y Af-body simulations (e.g., I White! [l986; 
: al.ll2006h : i 



Secondly, generic features of the e nsuing equilib- 
rium structures have been sought (see IHansenl 



Firstly, the halo growth has been recognized (see|Zhao et al. 


2003 


; Wechsler et al. 


2006; Hoffman et al. 2007; Diemand et al. 


2007 


) to comprise two stages: an early fast collapse including a 



few violent major mergers, that builds up the halo main 'body' 
with structure set by dynamical relaxation; and a later, quasi- 
equilibrium stage when the body is nearly unaffected, while 
the outskirts develop from the inside-out by min or mergers and 
smooth accretion (see Salvad or-Sole et al . 2007). The transition 
is provided by the time when a DM gravitational well attains 
its maximal depth, i.e., the radial peak of the circular velocity 
v 2 = G M/ R attains its m aximal height, along a given growth 
history (see lLiet al]l2007l) . 



; quilib- 
[20041: 



iDehnen & McLaughlin! 120051: ISchmidt et all 120081) among 
powerlaw correlations of the form cr 2 ^ 3 jp 2 ^ oc r a ; this involves 
the density p(r) and the velocity dispersi on cr 2 D = c 2 (1 + D j3), 
with anisotropy inserted via the standard Binnev ( 19781) param- 
eter (3 = 1 - (Tg/cr 2 and modulated by the index D (see IHansenl 
12007!) . It i s matter of debate w hich of these corre l ations best 
apply, see ISchmidt et ail (120081) and iNavarro et ail d2008l) : the 
former authors, in particular, find that the structure of different 
simulated halos may be described by different values of D, with 
linearly related values of e and a (see their Eqs. 4 and 5). 
Here we shall focus on the specific correlation 



K = 



P 



°7 

2/3 



(1) 



that involves solely the squared radial dispersion <x 2 , correspond- 
ing to D — and e — 3. This is because K has not only 
the striking form of a DM 'entropy' (or rather adiabat), but 
also the related operational advantage of providing a direct ex- 
pression of the radial pressure term per 2 — Kp 5 ^ oc r" p 5 ^ 3 
appearing in the Jeans equation for the radial equilibrium; in 
the latter any anisotropy is already accounted for by a sepa- 
rate term (see Eq. 2 below). On the other hand, the correla- 
tion K oc r" with a w 1.25-1.3 p rovides a simple yet ef- 
fective fit of ma n y simulations (see iTaylor & Navarro] [200 1 ; 
i Rasia et alj 12004 Hoffman et al J 120071; IDiemand et all 12 007; 
I Schmidt et alJl2008l;|Ascasibar & Gotfloberll2008!;lNavarro et al.l 
2008; IVass et al.ll2008l and many others). In the lower a range, 
Eq. (1) has the added bonus of preserving the classic self-similar 
slope in the halo body (see Eq. 3 below). 

To independently probe the matter, lLapi & Cavalierd 
(2009a) performed a semianalytical study of the two-stage halo 
development, and derived (consistently with the simulations) 
that a is set at the transition time via scale-free stratification of 
the particle orbits throughout the halo body, and thereafter re- 
mains closely constant and uniform at a value within the narrow 
range 1.25 - 1.3. Moreover, they found that on average the val- 
ues of a depend though weakly on the mass of the halo, such 
that a » 1 .3 applies to galaxy clusters, while a « 1 .25 applies to 
Milky Way sized galaxies. 



2.2. The DM a-profiles 

The halo physical profiles may be derived from the radial Jeans 
equation, with the radial p ressure per 2 , oc r" p 5 ^ and anisotropics 
described by the standard Binney (f 19781) parameter p. Thus the 
Jeans equation simply writes 



3 / v 2 \ 6 



(2) 



in terms of the l ogarithmic d e nsity s lope y = -dlogp/dlogr. As 
first sh own by lAustin et al.l ([2005) and IDehnen & McLaughlin 
(2005), Jeans supplemented with the mass definition M{< r) = 
4-n dr' r' 2 p(r') entering v 2 = GM(< r)jr, provides an integro- 
differential equation for p(r), that by double differentiation re- 
duces to a handy 2 nd -order differential equation for y. 

To set the context for the Milky Way DM distribution, we 
recall that the space of solutions for Eq. (2) spans the range a < 
1.296; the one for the upper bound and the behaviors of others 
one s have been analytically invest iga ted by lAustin et alTj[ 2005) 
and lDehnen &~ McLaughlin (2005). In lLapi & Cavalierd (12009a) 
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we explicitly derive the Jeans solutions with /? = (meaning 
isotropy) for the full range a r « 1.25 - 1.296 subjected to regular 
boundary conditions both at the center and in the outskirts, i.e., a 
round minimum of the potential with a finite pressure (or energy 
density), and a. finite (hence definite) overall mass, respectively. 
These we dubbed 'a-profiles' . 

The corresponding density runs steepen monotonically out- 
wards, and are summarized by the pivotal slopes 

3 3 

y = 6-3a, yb = ^(^+a); (3) 

these start from the central (r — > 0) value y a x 0.75 - 078, 
steepen in the halo main body to yo ~ 2.25 -2.1 (the former be- 
ing the slope from the classic self-similar collapse), and steepen 
further into the outskirts to typical values % ss 3.38 - 3.44 before 
a final cutoff. Thus the inner slope is considerably /faffer and the 
outer slope steeper c ompared to the empirical NFW formula (see 
iNavarro et al.lll997h ; in comparison to the Einasto profile, the 
main difference occurs in the inner regions where the ff-profile 
rather than fiat is moderately steep (see Eq. 3 and Fig. 1). 

For a density profile, a relevant parameter is the 'concentra- 
tion' c = R v /r-2, defined in terms of the virial radius R v and of 
the radius r_2 where y — 2; in the context of ff-profiles c may be 
viewed as a measure either of central condensation (small r_2) 
or of outskirts' extension (large R v ). The concentration consti- 
tutes an indicator of the halo age; in fact, numerical experiments 
(seelBullock et alJl200U IZhao et al.ll2003l: IWechsler etaT]|2006l 
iDiemand et al.ll2007l) show that c(z) ~ 3.5 holds at the end of the 
fast collapse stage, to grow as c(z) ~ 3.5 (1 + z,)/(l + z) during 
the slow accretion stage after the transition at z t - Current values 
c » 3.5 (1 + Zt) ~ 10 apply for a galaxy like the Milky Way that 
had its transition at Zt S> 2. 

The density and mass distribution in the Milky Way are il- 
lustrated in Fig. 1 for the isotropic a-profiles with a = 1.25 
(dashed), for the NFW formula (dotted), and for the Einasto pro- 
file (dot-dashed). All densities have been normalized to the lo- 
cal density 0.3 GeV cirT 3 at the Sun's location r w 8.5 kpc 
within the Galaxy. We further adopt r_2 = 20 kpc (consistent 
with c = 10). Note from Fig. 1 that the Einasto and NFW pro- 
files differ substantially at the center as for the density, and in the 
outskirts as for the mass, while the ff-profile strikes an interme- 
diate course between the two. 

2.3. Anisotropy 

It is clear from Eq. (2) that anisotropy will steepen the density 
run for positive /3 meaning radial velocity dominance, as ex- 
pected in the outskirts from infalling cold matter. On the other 
hand, tangential components (corresponding to /? < 0) must de- 
velop toward the center, as expected from increasing importance 
of angular momentu m effects. This view is supported by numer- 
ical simulations (seelAustin et al.ll2005l;lHansen & Moorell2006t 
iDehnen & McLaughlirj|l2005l) . which in detail suggest the effec- 
tive linear approximation 

P(r)«l3(Q)+P'[y(r)-ya] (4) 

to hol d with 0(0) > -0.1 and « 0.2, limited to f3(r) < 0.5. 

In lLapi & Cavalierel (|2009b) we extended the ff-profiles to 
such anisotropic conditions i n the full range a » 1 25 - 1.3, 
inspired by the analysis by IDehnen & McLaughlin! (120051) for 
the specific case a « 1.3. We find the corresponding p(r) to be 
somewhat flattened at the center by a weakly negative /3(Q), and 
further steepened into the outskirts where /3(r) grows substan- 
tially positive. Specifically, the following simple rules turn out 
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Fig. 1. Density and mass profiles in the Milky Way. The dashed 
and solid lines illustrate the ff-profiles with a = 1.25 in the 
isotropic and the anisotropic case, with y a = 0.75 and 0.63, re- 
spectively; the dotted line represents the NFW formula, and the 
dot-dashed line refers to the Einasto profile. All profiles are nor- 
malized to the local density 0.3 GeV cirT 3 at the Sun's location 
r Q * 8.5 kpc within the Galaxy; we have adopted r_2 = 20 kpc 
and c = 10, see § 2. 



to apply: the slope ft' i n Eq. (4) drops out from the de rivatives of 
the Jeans equation (see Dehnen & McLaughlin 2005); the upper 
bound to ff now reads 35/27 - 4/3(0)/27; moreover, y a is modi- 
fied into 3ff/5 + 6/5(0)/ 5 while yo and yi, retain their form. 

The anisotropic ff-profiles for the Milky Way are shown as 
solid lines in Fig. 1 . We note, in particular, that even a limited 
central anisotropy (e.g.,yS[0] * -0.1) causes an appreciable flat- 
tening of the inner density slope bringing it down to y a * 0.63 
for a; = 1.25. This, of course, results in an even more consider- 
able flattening for the slope of the squared density, the relevant 
quantity in our context of DM annihilations. 
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2.4. A guide to profile computations 



Finally, in the Appendix we provide user-friendly analytic fits 
for the density runs of the a-profiles in terms of standard de- 
projected Sersic formulas, but with parameters directly derived 
from the Jeans equation. 

We stress that these physical a-profiles with their analytic fits 
are relevant to, and recently tested in several contexts, includ- 
ing the interpretation of gravitational lensing observations (see 
lLapi & Ca valiere 2009b), the physics of the hot diffuse baryon s 
constituting the Intra-Clust er Plasma (see | Cayaliere et al.| [2009), 
and galaxy kinematics (see lLapi & Cavalierell2009cl) . In the fol- 
lowing we focus on the specific a-profile with a — 1 .25 suitable 
for the Milky Way halo (see § 2. 1) to predict the DM annihilation 
signal from the GC. 



3. y-ray signal from DM annihilation at the GC 

The y-ray flux per solid angle due to DM annihilation along a 
direction at an angle if/ relative to the l.o.s. toward the GC may 
be written (under the commonly assumed spherical symmetry) 
as 



d0> 7 
dll 



= 3.74 x 10" 6 ^ 



<Iv> 



io- 



50GeV 



) 2 m (5) 



in units of m~ 2 s sr . The above expression is naturally 
factorized into a microscopic and an astrophysical term (e.g., 
Bergstrom 2009, and references therein). The former involves 
the mass of the DM particle wdm, the number of photons N 7 
created per annihilation, and the angle-velocity averaged anni- 
hilation rate (Sv) in terms of the particles' cross section S and 
velocity v. 

For the sake of definiteness we begin from considering a 
neutralino DM particle with mass Mdm ~ 50 GeV, annihilat- 
ing through the bb channel (with 100% branching ratio). We use 
the benchmark value for the annihilation rate (Sv) ~ 3 x 10~ 26 
cm 3 s _1 , corresponding to a thermal relic with a density close to 
the cosmological DM abundance 



3 x 10- 



<Iv> 



0.1 



(6) 



as measured by WMAP (see lDunkley et al.1l2009l) . To compute 
N y — J dE dNy/dEy we adopt a photon annihilation spectrum 
with shape 



rjx e 



b+cx+dx 2 +e.x 3 



(1) 



dN Y 
dx 

obtained from extrapolating the results bv lFornengo et al.l f2004) 
down to energies E y as 200 MeV; here x = Eyfm^M is the energy 
normalized to the DM mass, while rj = 1, a = -1.5, b = 0.579, 
c = -17.6080, d = 23.862, e = -25.181 are fitting parameters 
for the adopted microscopic DM model (see above). 

The astrophysical term of Eq. (5) is given by the integral of 
the (squared) DM density projected along the l.o.s. 



J r e 



dl p 2 (r) 



r© p 2 (r Q ) ' 



(8) 



normalized to p(r Q ) ai 0.3 GeV cm 4 , the local density at the 
Sun's location r Q ai 8.5 kpc. A non-trivial angular dependence 
results from the peripheral position of the Sun within the Milky 
Way halo, and involves only the angle ifi between the observed 
direction of the sky and the GC; in terms of Galactic latitude 



o 

cn 



cn 
o 



DM Annichik 


tion Signa 








1 




N FW 


in 


\ 


- a=1,25 iso 


cn 


'a 


— a= 1 .25 aniso " 
Einasto 


-1 7 






-2 a 
~o 






\ 
^ e " 

cn 






o 






-4 



50 100 150 

1> [ de g] 




: 1.0 




: 0.5 




: 0.0 


cn 


: -0.5 




: -1.0 


a 


; -1.5 




: -2.0 


CP 

o 



0.1 



1.0 
ifj [deg] 



10.0 



Fig. 2. The astrophysical factor J (if/) normalized to the value at 
if/ — 90 deg (left axis), and the corresponding annihilation flux 
per unit solid angle for E y > 200 MeV (right axis; see § 3 for 
details). The bottom panel zooms on the inner angular scales. 



b and longitude / with cos if/ = cos b cos /, the radial variable 
can be expressed as r — (r| + i 2 — 2 r { cos ifi) 1 ^ 1 on using the 
distance I along the l.o.s. Finally, when observing a region at an 
angular resolution AQ, one has to consider the average value of 
J, namely, 



f dQ J{iff) 
Jaci 



(9) 



with d£2 = cos b db dl. 

We compute and report in Table 1 the values of J at angu- 
lar resolutions AQ. = 10~ 3 sr and 10~ 5 sr for the a-profile with 
a — 1.25 in the isotropic and anisotropic cases, for the NFW 
formula, and for the Einasto profile. In Table 2 we list the cor- 
responding values of the y-ray flux for energies E y > 200 MeV. 
These outcomes are illustrated in Fig. 2. It is seen that relative 
to the NFW distribution, the fluxes predicted from the isotropic 
and anisotropic a-profile are lower by factors from a few to sev- 
eral. Such fluxes are still within the reach of the Fer mi satellite; 
in fact , on the basis of the s imulations performe d bv lBaltz et al.l 
( 20081 IStrianl (l2009h . and IVitale et ail d2009al) . we expect the 
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Table 1. Values of the astrophysical factor J. 



DM profile 


An = i(r J 


An = io -5 


NFW 


1.2 x 10 3 


1.05 x 10 4 


a iso 


1.3 x 10 3 


6.5 x 10 3 


a aniso 


5.1 x 10 2 


1.6 x 10 3 


Einasto 


1.1 x 10 2 


1.3 x 10 2 



Table 2. Values of the y-ray flux (in rrT 2 s" 1 ) for E y > 200 MeV. 



DM profile 


An = 10 


-3 


An = io- 5 


NFW 


4.7 x 10 


-4 


4.0 x 10" s 


a iso 


5.2 x 10 


4 


2.5 x 10" s 


a aniso 


2.0 x 10 


4 


6.0 x 10- 6 


Einasto 


4.1 x 10- 


-5 


5.0 x 10~ 7 



annihilation signal to be probed at a 3-<x confidence level over a 
few years. 

The above values may be compared with the current upper 
bound to the integrated flux of 2.43*^ x 10~ 3 irT 2 s -1 based 
on Fermi measurements at E y > 200 MeV during 8-month ob- 
servations of the GC over a solid angle AQ » 10~ 3 sr (see 
lAbdo etal.ll2009t lAtwood et al.ll2009t Ivitale etalJl2009bh : this 
bound decreases as 1 / '\fi with the observation time t . However, 
the flux currently observed includes contributions from diffuse 
or not yet resolve d Galactic sources, that are being progres- 
sively removed (seelGoodenough & Hooper 2009; Striani 2009; 
Vital e et alj l2009bl) ; next stages of such a process will take 
longer observations aimed at determining the spectrum of indi- 
vidual resolved sources and a c areful likelyhood ana lysis of the 
backgrounds (see discussion bv lCesarini et al.ll2004j) . 



4. Discussion and conclusions 

We have presented our a-profile with a — 1 .25 for the equilib- 
rium density and mass distributions in a galactic DM halo, and 
specifically in the Milky Way. We have shown that this profile 
constitutes the robust solution of the equilibrium Jeans equation 
with physical inner and outer boundary conditions, i.e., finite 
pressure and round potential minimum at the center, and finite 
overall mass. The corresponding density profile p(r) is intrin- 
sically flatter at the center, and intrinsically steeper in the out- 
skirts, relative to the empirical NFW formula. These features are 
sharpened yet in halos with anisotropic random velocities. We 
have also provided the reader with a precise and user-friendly 
analytic fit to the a-profile (see Appendix for details). 

Then we have focused on the role of this a-profile as a bench- 
mark for computing the DM annihilation signal expected from 
the GC. In fact, we have computed the 'astrophysical factor' 
J(i[f) (angular distribution, independent of microphysics) enter- 
ing the expression of the annihilation flux. As a definite exam- 
ple, we have also computed the y-ray flux on adopting a simple, 
fiducial microscopic model; this we find consistent with current 



Fermi observations, given that the latter may include contribu- 
tions from still unresolved point sources. 

Given the physical a-profile and the corresponding factor 
J(if/), the extension to m ore complex microscopic scenarios like 
mSUGRA (started by [Chamseddine et all 119821: iBarbieri et all 
119821; lOhta et all 119831; lHall et all 119831) will be easily made 
in terms of annihilations channels, cross sections and particle 
masses. In this context our a-profile relieves astrophysical un- 
certainties related to the macroscopic DM distribution. We stress 
that constraints on particle cross sections and masses inferred 
from radio and y-ray observations of the GC have been to now 
more sensitive to the assumed DM distribution than to specific 
annihilati on channels ( d ifferen t from leptonic rf), see Figs. 3 
and 4 in Bert one et al.l d2009l) . In fact, the latter show that a 
DM distribution with an inner slope like our a-profile is re- 
quired to allow cross sections (S v) ^ 10~ 26 cm 3 s _1 with masses 
;«dm ~ 500 GeV for the non-leptonic channels that are widely 
considered on grounds of theoretical microphysics. 

Concerning small scales r ~ a few tens of pes around the GC, 
we touch upon a number of possible deviations of the very in- 
ner DM density distribution from our benchmark a-profile (solid 
line in Fig. 1, top panel). For example, the process of galaxy for- 
mation could lead either to flattening or to some steepening of 
the inner DM distribution. The former may occur either owing to 
transfer of e nergy and/or angular momentum from th e baryons to 
the DM (see lEl-Zant et alj200UlTonini et al]|2006l) . or owing to 
quick mass removal following the energy fe edback from stars 
or active galactic nuc lei (see discussion by lLauer et al J 12007: 
iKormendv et al.ll2009l) . On the other hand, steepening might be 
induced by the 'adiabatic' contraction of the b aryons into the 
disc (see lBlumenthal et al.ll986tlMo et aljl998b ; but even in ex- 
treme cases (see discussion bv lAbadi et al.l20 09) such a contrac- 



tion would yield an inner DM density profile p(r) oc r 



-3/(4-y„) 

still flatter than 1 though somewhat steeper than the original 
y a ~ 0.63 - 0.75. Fin ally, at the very cent er of the Gala xy any 
accretion of DM (e.g jGondolo et all 1 999t iBertone et al.f 2002) 
onto the nuclear supermassive black hole might enhance the DM 
distribution on tiny scales r < 10 _I pc. 

Summing up, we stress that all such alterations of the inner 
slope would occur on scales smaller than some 10 pes; although 
significant at levels of a few percent to account for the central 
stellar lighQ, their import is far smaller for what the annihilation 
signal is concerned, and in the average over 10 _1 deg the flux is 
altered by less than 0.1%. In fact, these corrections are currently 
at, or below the resolution limit and the prospective sensitivity 
of Fermi. 

Other possible targets include the dwarf spheroidal galaxies 
in the Local Group. These on the one hand constitute cleaner en- 
vironments than the GC owing to their dearth of stellar sources; 
on the other hand, their distance if modest on intergalactic scales, 
already makes detecting and resolving the related anni hilation 
signal a real challenge for Fermi (e.g.. |Pieri et al .1120091) . In ad- 
dition, the shallow gravitational potential wells of these systems 
make them particularly prone to energy feedback events (see 
above), that may flatten the inner DM distribution to flat slopes 
y a < 0.63 (consistent with kinematical observations), to the ef- 
fect of further lowering the annihilation signals. Upper limits 
more stringent than the current value (Ev) < 10~ 25 cm 3 s _1 at 
a mass mdm ~ 50 GeV will require delicate stacking over an 
ensemble of dwarfs. 



1 We note that a flat slope is known to describe the very central light 
distribution in luminous ellip t icals, related to compl ex small-scale dy- 
namics (see lLauer et al.ll2007l : lKormendv et al]|2009l) . 
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Concerning particle cross section and masses, we recall 
that the PAMELA satellite recently observed an excess of the 
positron fraction e + /(e + + e~) in the cosmic ray e ± spectra rel- 
ative Jo_Jhe_e2cp_ected_astrophysical background above 10 GeV 
(see lAdriani et al. 2009). This excess can be simply explained 
in terms of a single or a few sources like pulsars, that are ex- 
pected to produce a powerlaw spe ctrum of e ± pairs with a cutoff 
at several TeVs (see lBertonel20 09) . On the other hand, the signal 
may be also interpreted in terms o f DM annihilations occurring 
throughout the Galactic halo (e.g jBertone et al. 2009). 

If this is to be the case, however, the flux measured by 
PAMELA mandates for very large effective annihilation cross 
sections (Ev) ~ 10~ 23 cm 3 s _1 , well above the natural value 
suggested by the cosmological DM abundance (see § 3). From 
a microphysical point of view, this is still conceivable in 
scenarios with Sommerfe ld enhancements (see discussion by 
lArk ani-Ha med et al. 2009); the cross section may be enhanced 
by a factor ~ 10 2 for velocities v/c ~ 10~ 3 . On the other hand, 
such a large Sommerfeld effects would also yield a strong y-ray 
annihilation signal towards the GC; for this, little room is al- 
lowed on the basis of the current upper limit provided by Fermi 
(see § 3), unless the DM particle mass substantially exceeds 50 
GeV. 

Another possibility is to invoke a large boost factor of the ef- 
fective cross section due to dumpiness in the Galactic halo, i.e., 
a crowd of dense subhalos; however, state-of-the-art numerical 
simulations suggest such boosts not to be rea listic in the Galaxy, 
even less at the GC ( see lspringel et aT1l2008l and discussion by 
lLattanzi & Siikll2009l) . 

To sum up, we have discussed why the a-profile with a = 
1.25 (see Fig. 1) constitutes a reliable DM distribution in the 
Galaxy; we have argued that it will provide a benchmark to 
gauge in terms of DM annihilation the y rays from the GC to 
be detected with Fermi (see Fig. 2). Such an a-profile will be 
instrumental to derive reliable information concerning the mi- 
croscopic nature of the DM particles. 
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Appendix A: Analytic fit to the a-prof ile 

To complement t h e analy tical details extensively dealt with by 
lLapi & Cavaii ere (2009 a|) and to enable a straightforward com- 
parison with the classic NFW and Einasto density runs, here we 
provide a handy analytic fit to the a-profiles in terms of the de- 
projected Sersic formula substantiated with parameters directly 
derived from the Jeans eq uation. We base on the expression (see 
iPrugniel & Simienll 19971) 



Table A.l. Values of the fitting parameters of Eq. (Al) in the 
isotropic case; a — 1 .25 applies for the Galactic halo. 



^Uf-I cxp<- 
P(r-i) I r- 2 ' ' ] 



2-t 



r 

r-2 



- 1 



(A.l) 



where t and rj are two fitting parameters; the standard Einasto 
profile obtains for t = 0. The values of t and r\ for different a 
of interest here are reported in Tables A.l and A. 2 both in the 
isotropic and the anisotropic cases; note that t m y a is required 
by the physical boundary condition satisfied at the center (see 
§ 2). The resulting fits to the density runs of the ar-profiles hold 
to better than 10% in the relevant range 



a 1.25 1.26 1.27 
r 0.750 0.756 0.762 
j) 0.319 0.298 0.277 



Table A.2. Values of the fitting parameters of Eq. (Al) in the 
anisotropic case; a — 1.25 applies for the Galactic halo. 



a 1.25 1.26 1.27 
r 0.630 0.636 0.642 
ri 0.364 0.342 0.319 



The mass corresponding to the density distribution of 
Eq. (Al) reads 



M(< r) 



3-t 2-t / r 
T] ' T] \r_ 2 



(A.2) 



where T[a,x] = f*dt f~ x e"'/j£* dtf 1 e 'is the (normalized) 
incomplete F-function. 
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